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1.  Introduction 


Constraint  Logic  Programming  languages  span  the  void  between  symbolic  and  minieric 
computation.  The  increased  power  of  expression  that  these  languages  present  to  a  pro¬ 
grammer  expands  the  field  of  applications  that  logic-based,  declarative  languages  can 
represent.  Tliis  report  is  an  initial  effort  to  introduce  the  specific  language  CLP(iK)  to 
the  planning  and  scheduling  area  of  artificial  intelligence. 

Section  2  of  this  work  is  an  overview  of  Constraint  Logic  Programming,  including  the 
motivation  and  design  goals  of  CLP(9t).  Short  examples  are  also  given  in  this  section 
to  illustrate  new  and  useful  features  of  the  language.  A  description  of  the  CLP(5i) 
architecture  is  added  to  complete  the  system’s  basic  structure.  Finally,  this  section 
presents  three  application  areas  of  CLP(9?):  (1)  financial  analysis.  (2)  problems  from 
electrical  engineering,  and  (3)  analysis  and  partial  synthesis  of  truss  structures. 

Section  3  describes  two  current  planning  or  scheduling  systems.  The  short  summaries 
that  are  presented  there  are  intended  to  emphasize  the  basic  features  or  to  illustrate  those 
characteristics  of  the  systems  that  are  deemed  appropriate  for  CLP(9\). 

Section  4  discusses  two  original  planning  examples  that  link  CLP(9?)  with  planning 
problems.  The  codes  for  the  two  examples  are  contained  in  the  appendices.  Conclu¬ 
sions  for  this  work  are  found  in  Section  5. 


2.  Constraint  Logic  Programming 


"The  class  of  Constraint  Logic  Programming  tCLP)  languages  has  caused  much  ex¬ 
citement  in  computer  science  because  of  its  conglomeration  of  such  diverse  areas  as 
numerical  analysis,  artificial  intelligence,  operations  research,  logic,  mathematics,  and 
formal  languages"  [Cohen  1990].  The  combination  of  logic  programming  with  con¬ 
straint  satisfaction  is  a  natural  extension  of  traditional  logic  programming  languages 
since  the  fundamental  unification  process  can  be  viewed  as  a  special  case  of  constraint 
satisfaction;  i.e.,  unification  of  f(X,a)  with  f(b,Y)  is  equivalent  to  solving  f(X.a)=f(b,Y). 

Constraint  Logic  Programming  (CLP)  is  a  class  of  logic  programming  languages  that 
are  modeled  on  rule-based  constraint  programming.  The  specific  instant  CLPCiH)  uses 
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the  real  numbers  as  its  structure  of  computation.  In  panicular,  ihc  consirainis  arc 
relations  over  "real  arithmetic  terms  with  umnterpreted  functors'  [Jaffar  8o].  Since  all 
CLP  languages  are  soundly  based  on  a  single  framework  of  formal  semantics,  there  are 
no  restrictions  in  CLPCtK)  to  the  Herbraiid  I  niverse,  to  unification,  or  to  allowed 
equations. 

There  are  several  features  of  CLP(9?)  that  are  note-worthy.  The  most  important  features 
are  briefly  described  here  and  illustrated  more  fully  later.  The  output  of  CLPtrK)  pro¬ 
grams  is  symbolic;  thus,  it  is  not  only  possible  to  see  responses  to  queries  like  X=o  or 
Y=bill,  but  one  may  also  have  X*X  =  Y+d. 

CLP(9l)  also  expands  the  use  of  declarative  programming,  for  example,  in  using  an 
equation  of  the  form  X  =  Y  4-  4,  it  is  no  longer  necessary  to  wnie  two  different  rules, 
depending  on  whether  X  or  Y  is  instantiated  at  the  time  that  the  equation  is  invoked. 
This  equation  is  merely  a  constraint,  and  CLP(9?)  is  able  to  deal  wuh  all  four  possible 
instantiation  cases:  neither  X  nor  Y  is  ground.  X  is  ground.  Y  is  ground,  or  both  and 
Y  are  ground.  From  the  programmer’s  view,  it  is  only  required  to  state  the  equation. 
CLP(S^)  then  makes  the  appropriate  interpretation  during  execution. 

CLP(S^)  is  able  to  run  Prolog  programs  with  no  changes.  Indeed,  one  of  the  design 
goals  of  CLPCiR)  was  not  only  to  execute  Prolog  programs  but  to  have  them  run  as 
efficiently  as  normal.  This  required  separating  the  interpreter  into  an  inference  engine 
and  a  constraint  solver.  In  running  a  Prolog  program,  there  should  be  no  need  for  the 
constraint  solver.  Consequently,  the  inference  engine  is  essentially  a  Prolog  interpreter. 

The  second  fundamental  design  goal  of  CLP(5\)  dictated  that  only  linear  constraints 
would  be  considered  by  the  constraint  solver.  This  limitation  was  ameli^oated  by  de¬ 
fining  a  delay  and  wakeup  condition.  Basically,  non-linear  constraints  p;e  delayed  until 
a  sufficient  number  of  the  variables  are  instantiated  and  to  make  the  constraints  linear. 

The  third  design  goal  required  CLP(9?)  to  deal  with  constraints  incrementally.  This 
means  that  in  the  execution  of  a  query  it  is  not  necessary  to  v'erify  the  solvability  of  an 
entire  set  of  constraints  when  a  single  new  constraint  is  added.  Rather  it  is  necessary 
only  to  verify  that  the  new  constraint  is  compatible  with  the  old  set  of  constraints. 

Other  design  goals  of  CLP(5l)  included  a  choice  of  appropriate  algorithms  for  testing 
the  satisficability  of  systems  of  constraints,  a  demand  that  the  system  be  able  to  simplify 
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a  system  of  constraints,  and  that  the  constraints  be  represented  internally  in  a  canonical 
form.  These  latter  objectives  deal  with  the  very  practical  issues  of  constructing  an 
efficient  implementation.  They  will  not  be  examined  funher,  but  their  solution  will  be 
lauded. 


2.1  Constraints  versus  Assignments 

In  this  section  two  examples  are  given  to  illustrate  the  increased  power  of  expression 
available  to  CLP(91).  Both  examples  involve  replacement  of  the  assignment  statement 
("is")  in  Prolog  with  a  constraint. 

Consider  the  following  predicate  which  defines  the  addition  of  the  first  two  argument: 

plus(X.Y.Z) Z  is  X+Y. 

This  rule  assumes  that  the  variables  X  and  Y  have  been  instantiated  at  the  time  the  rule 
is  invoked.  If  either  X  or  Y  is  not  ground  or  if  Z  is  ground  at  invocation,  a  run-tmie 
error  occurs.  In  order  to  make  this  rule  valid  in  all  possible  circumstances,  it  is  neces¬ 
sary  to  replace  it  with  the  following  set  of  rules: 

plus(X,Y,Z) not(var(X)),  not(var(Y)),  not(var{Z)), 

S  is  X+Y,  Z=S. 

plus(X,Y,Z) not(var(X)),  not(var(Y)),  var(Z).  Z  is  X+Y. 
plus(X,Y,Z) not(var(X)),  var(Y).  noi(var(Z)),  Y  is  Z-X. 
plus(X,Y.Z) var(X),  not(var(Y)),  not(var(Z)),  X  is  Z-Y. 

CLP(9l)  allows  this  set  of  rules  to  be  replaced  by  this  more  declarative  representation: 

plus(X,Y,Z) Z=X+Y. 

Here  the  intention  is  much  clearer;  the  variables  X.  Y,  and  Z  are  constricted  to  obey  the 
given  equation.  It  is  not  necessary  for  the  programmer  to  described  separately  the  four 
instantiation  cases. 

Moreover,  it  is  equally  awkward  in  traditional  Prolog  to  consider  the  other  four  possi¬ 
bilities  of  variable  instantiations:  two  of  the  three  variables  are  uninstantiated  or  all 
three  variables  are  uninstantiated.  CLP(9^)  accepts  the  single  equation  and  internally 
considers  the  separate  cases.  In  CLP(SR)  the  response  may  be  either  yes/no  (when  all 
variables  are  instantiated),  an  assignment  (when  two  of  the  three  variables  are  instanti¬ 
ated),  or  an  equation  (in  other  cases). 


The  second  example  compares  the  definition  of  the  Fibonacci  senes  in  Prolog  and  in 
CLPC31). 


in  Prolog 
fib(OJ). 
fib(i,l). 
fib(N.R) 

N1  is  N-1, 


m  CLPm) 
fib(0,l). 
fib{l,i). 
fib(N,Rl+R2) 
N>=2, 


fib(Nl,Rl), 
N2  is  N-2, 
fib(N2,Rl), 


fib(N-l,Rl). 
fib(N-2,R2). 
R  isRl4-R2. 


Both  of  these  programs  can  answer  the  question  ?-fib(  10. F)  with  the  answer  F~89.  but 
only  the  CLP(SR)  program  can  answer  the  question  ?-fib(N.89}  with  N-10.  Tliis  ex¬ 
ample  illustrates  several  features  of  CLP(9l).  First,  the  second  program  is  more 
declarative:  to  find  the  Nth  Fibonacci  number,  add  the  Fibonacci  numbers  of  the  (N-1  ).st 
and  (N-2)nd  Fibonacci  numbers  if  N  is  larger  than  one.  Secondly,  the  CLP(9\)  defini¬ 
tion  allows  a  reversal  of  role  of  the  "input"  variable  and  the  "output"  variable.  Although 
this  is  common  with  many  other  Prolog  predicates  (compare,  the  standard  definitions  of 
member  and  append),  it  is  not  the  intention  here.  Finally,  the  CLP(9^)  version  allows 
arithmetic  operations  to  occur  in  the  arguments  of  a  predicate.  This  is  a  radical  notion 
only  for  Prolog  since  most  procedural  languages  have  generally  provided  such  facilities 
for  procedure  calls. 


2.2  The  CLP(5R)  Interpreter 

Arithmetic  terms  in  CLP(9I)  are  constructed  from  real  constants  .variables,  and  these 
interpreted  functors:  +,  *  /,  sin,  cos,  tan,  pow  (the  power  function),  min,  and  max. 

The  CLP(9t)  interpreter  is  composed  of  five  pans:  (1)  the  inference  engine.  (2)  the  in¬ 
terface  between  the  inference  engine  and  the  constraint  solver,  (3)  the  equation  solver. 
(4)  the  inequality  solver,  and  (5)  the  output  module. 

The  constraint  solver  itself  is  composed  of  the  equation  and  the  inequality  solvers.  The 
inference  engine  is  basically  a  Prolog  interpreter,  which  controls  the  execution  of  the 
derivation  steps  and  updates  variable  bindings.  The  usual  stack  mechanisms  and  sym¬ 
bol  tables  are  maintained  here.  The  interface  module's  primary  purpose  is  to  smooth 
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the  transition  between  the  inference  engine  and  the  pieces  of  the  consn-aini  solver.  The 
interface  evaluates  complex  arithmetic  expressions  and  transforms  constraints  into 
standard  forms. 

Constraints  that  are  received  by  the  constraint  solver  occur  in  one  of  three  forms;  linear 
equations,  linear  inequalities,  or  non-linear  equations.  Linear  equations  are  sent  to  the 
equation  solver  while  linear  inequalities  are  sent  to  the  inequality  solver.  Non-lmear 
equations  are  delayed  until  so  many  of  its  variables  are  ground  that  the  equations  be¬ 
come  linear.  The  output  module  finally  converts  the  internal  representations  of 
constraints  to  a  canonical  output  form. 

"The  central  data  structure  in  the  equality  solver  is  a  tableau  which  stores,  in  each  row. 
a  representation  of  seme  variables  in  the  form  of  a  linear  combination  of  parametric 
variables"  [Jaffar  1988]. 

As  progress  is  made  with  the  linear  constraints,  the  involved  variables  become  ground. 
It  is  then  possible  to  check  if  these  variables  are  involved  in  any  of  the  delayed  non¬ 
linear  constraints.  After  making  appropriate  adjustments  to  the  non-linear  constraints, 
a  queue  is  formed  of  all  non-linear  constraints  that  have  become  linear.  These  linear 
constraints  are  then  passed  to  either  the  equation  solver  or  the  inequality  solver  before 
proceeding  with  the  next  derivation  step  of  the  proof.  Chaining  effects  in  which  one 
linearized  constraint  causes  another  non-linear  constraint  to  become  linear  are  possible. 

A  simple  bookkeeping  device  is  used  to  determine  when  backtracking  in  the  constraint 
solver  is  necessary.  In  general,  backtracking  in  the  constraint  solver  always  corre¬ 
sponds  to  a  backtracking  point  in  the  inference  engine  but  not  vice  versa. 

Linear  equations  are  stored  in  parametric  form.  Thus,  a  constraint  like  X  =  Y.  which 
involves  program  variables  X  and  Y  may  be  represented  as  these  two  parametric  equa¬ 
tions;  X  =  t  and  Y  =  t.  In  general,  any  variable  will  be  converted  to  this  parametric  form: 
X  =  b  -f  Citi  -I-  C2t2  +  +  Cictk. 


Linear  inequalities  are  solved  by  the  Two-Phase  Simplex  algorithm.  Details  of  this 
algorithm  are  given  in  [Jaffar  1986]. 
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2.3  Applications  of  CLP(5K) 


CLP(9^)  has  been  used  for  an  ever  increasing  number  of  appiicaiions.  These  three  ex¬ 
amples  are  cited  to  display  the  diversity:  (1)  options  trading  on  the  stock  market.  (2) 
electrical  engineering  circuit  analysis,  and  (3)  truss  analysis  and  panial  synthesis  ^)f 
structures  in  civil  engineering. 

2.3.1  CLP(9t)  in  Financial  Analysis 

Options  are  used  to  express  confidence  levels  in  the  intrinsic  value  of  an  asset.  For 
example,  one  may  make  a  call  option  to  buy  a  stock  before  some  specific  date.  Should 
the  stock  drop  very  low  before  the  specified  date,  the  owner  of  the  call  option  would 
typically  purchase  the  stock  at  its  low  point.  In  the  worst  case,  the  stock  increases 
monotonically,  and  the  owner  must  buy  the  stock  when  the  specified  date  arrives. 

Two  features  of  options  trading  have  been  cited  [Lassez  1987]  that  make  it  an  appro¬ 
priate  application  of  combining  symbolic  and  numeric  computation.  They  are  that 
traders  typically  combine  (1)  their  personal,  expert  heuristics  with  (2)  numeric  valua¬ 
tion  functions  in  decision-making. 

The  indifference  of  rules  in  CLP(9^)  to  dictating  which  variables  are  input  and  which 
variables  are  output  is  an  attractive  characteristic  to  option  trading  since  it  is  not  nec¬ 
essary  to  write  additional  rules  for  reversing  the  direction  of  computation.  Funhermore. 
in  CLP(9l)  it  is  possible  to  unite  symbolic  output  and  constraints  on  goals  to  obtain 
"what  if’  analyses. 

An  interesting  feature  of  option  trading  is  that  payoffs  for  option  combinations  can  be 
represented  by  piecewise  linear  functions.  These  functions  can  be  expressed  in  term  of 
two  elementary  functions:  (1)  the  Heaviside  function,  and  (2)  the  ramp  function,  as 
shown  in  [Lassez  1986].  The  Heaviside  function,  h(X,Y),  is  zero  if  Y  <  X  and  one  if  Y 
>=  X.  The  ramp  function,  r(X,Y),  is  zero  if  Y  <  X,  and  Y-X  if  Y  >=  X. 

2.3.2  CLP(St)  in  Electrical  Engineering 

A  second  area  of  applications  using  CLP(S^)  has  been  with  some  electrical  engineering 
problems  [Hemtze  1986].  Three  categories  of  these  problems  are  lumped  circuits,  dig- 
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ital  filtering  networks,  and  electromagnetic  fields. 

Circuit  problems  include  both  analysis  and  synthesis  of  analogue  circuits,  in  which  the 
constraints  are  either  local  (e.g..  Ohm’s  law  for  a  resistor)  or  global  (e.g.,  Kirchoff’s  law' 
in  nodal  analysis).  Specific  examples  of  these  problems  are  analysis  of  steady-state 
RLC  circuits  and  diode  components  in  RLC  circuits.  Bipolar  junction  transistors  for 
amplifier  and  logic  circuits  have  also  been  analyzed  and  designed. 

In  the  area  of  digital  signal  flow,  the  simulation  of  linear  shift-invariant  digital  system 
with  the  independent  variable  time  have  been  successfully  modeled  in  terms  of  linear 
signal-flow  graphs.  In  these  graphs  signal  values  at  the  nodes  at  successive  time  steps 
are  simulated  by  summing  constraints  for  each  node  and  by  collecting  branch  equations. 
In  CLP(9l)  one  needs  only  to  declaratively  state  the  equations  while  the  interpreter  de¬ 
cides  the  particular  method  for  solving  the  equations. 

The  final  category  of  electrical  engineering  problems  involves  numerical,  approximate 
solutions  for  the  pertinent  partial  differential  equations  of  electromagnetic  fields.  .■X 
finite  difference  approximation  is  the  common  approach  for  these  problems.  The  ex¬ 
amples  of  this  category  illustrate  the  elegance  of  CLP(iK)  in  using  finite  difference 
methods  by  its  concise  representations  and  its  use  of  the  constraint  collection 
mechanism. 

2.3.3  CLP(9^)  in  Truss  Structures 

Analysis  and  synthesis  problems  are  common  to  all  areas  of  engineering.  Recently 
[Lakmazaheri  1989]  has  shown  that  constraint  logic  programming  is  suitable  for  the 
analysis  and  partial  synthesis  of  truss  structures.  Typically  these  two  problems  of  anal¬ 
ysis  and  synthesis  are  considered  antithetical  in  engineering;  however,  a  CLP(iK) 
representation  is  able  to  use  only  one  representation  for  the  two.  The  particular  query' 
dictates  the  appropriate  form. 

This  work  on  truss  structures  has  shown  that  the  technique  of  local  propagation  is  in¬ 
adequate  for  satisfying  simultaneous  constraints  when  the  associated  constraint  graph 
contains  cycles.  But  in  the  constraint-based  approach  the  constraints  can  be  viewed  as 
relations  rather  than  functions;  hence,  there  is  no  explicit  commitment  to  a  specific 
variable  being  either  input  or  output. 


The  constraints  associated  with  structural  components  are  classified  by  these  cateeo- 

ries:  (1)  truss  components,  which  satisfy  a  two  or  three  dimensional  form  of  the  vector 

equation  KD=F  with  K  an  n  by  n  (n=2  or  3)  stiffness  matrix,  D  the  n  by  1  displacement 

vector  at  the  two  end  nodes  and  F  the  n  bv  1  force  vector  at  the  two  end  nodes;  (2) 

¥ 

support  components  of  the  type  pin  support,  rollerX  support,  or  lollerY  support;  and  (3) 
load  components,  each  of  which  is  associated  with  a  force  vector  and  a  displacement 
vector. 

The  position  and  displacement  vectors  are  functions,  and  the  nodes'  force  histories  are 
lists.  Nodes  in  two  dimensions  are  terms  of  this  form; 

n(p(Px,Py),  d(Dx,Dy),  H) 

where  Px  and  Py  are  the  x  and  y  coordinates  of  the  position  vector,  Dx  and  Dy  are  the 
X  and  y  coordinates  of  the  displacement  vector,  and  H  is  the  node  force  history  list. 

One  can  then  describe  general  rules  for  the  local  nodal  behaviors.  This  is  analogous  to 
Ohm’s  law  in  electrical  circuits.  For  example,  this  rule 

A  valid  node  exists  if  its  coordinates  are  valid. 
is  represented  by  this  CLP(SK)  rule 

node(n(p(Px,Py),_,_)) ;-  ordinate(Px),  ordinate(Py). 

Since  the  behavior  of  a  structural  model  is  derived  from  the  behavior  of  its  structural 
components,  rules  can  be  added  that  insure  the  satisficability  of  the  constraints  associ¬ 
ated  with  the  structural  components. 

Four  examples  are  presented  in  [Lakmazaheri  1989].  The  first  example  illustrates  the 
capability  of  the  constraint  logic  program  to  perform  numerical  analysis  of  a  statically 
indeterminate  truss  strucmre.  The  second  example  shows  the  capability  of  the  program 
to  perform  symbolic  computation  for  structural  analysis.  The  third  example  gives  a 
partial  synthesis  of  a  simple  truss  structure  using  the  same  program.  Finally,  the  last 
example  demonstrates  the  applicability  of  the  program  for  solving  more  complex 
problems. 

The  main  contribution  of  this  work  is  the  development  of  a  general  and  declarative 
constraint-based  formulation  for  the  analysis  and  partial  synthesis  of  truss  structures 
using  the  framework  of  constraint  logic  programming  methodology.  A  single  and  uni¬ 
form  framework  for  modeling  engineering  design  problems  facilitates  the  integration  of 
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different  design  activities  such  as  synthesis,  analysis,  and  evaluation. 

Currently  [Lakmazaheri  1992]  is  examining  two  strategies  to  improve  the  processing 
time  of  constraint-based  engineering  problems  by  using  distributed  processing  through 
problem  decomposition.  The  set  of  connectivity  constraints  then  ensures  displacement 
compatibility  and  force  equilibrium  on  the  connections  between  each  pair  of 
substructures.  A  study  of  5000  elements  has  shown  a  speed-up  factor  of  6  is  achieved 
when  the  structure  is  decomposed  into  sixteen  substructures. 


3.  Planning 

Constraints  are  a  familiar  part  of  planning  and  scheduling  problems.  The  choices  tha: 
are  made  in  such  a  design  involve  relations  (constraints)  among  objects  of  the  domain 
of  discourse.  The  use  of  constraints  in  anificial  intelligence  was  proposed  in  [Steele 
1980]  and  [Sussman  1980], 

This  section  describes  briefly  two  current  planning  or  scheduling  systems.  The  short 
summaries  that  are  presented  here  are  intended  to  emphasize  the  basic  features  or  to 
illustrate  those  characteristics  of  tire  systems  that  are  deemed  appropriate  for  CLP(S^). 

3.1  OPIS 

The  difficulties  in  factory  scheduling,  as  defined  in  [Smith  1990]  are  caused  by  the 
complexity  of  assigning  shared  resources  to  competing  processes  and  by  the  unpre¬ 
dictability  of  maintaining  a  stable  schedule  of  good  quality  within  acceptable  time 
limits,  opportunistic  Intelligent  Scheduling  (OPIS)  is  an  incremental  scheduler  that 
attacks  these  difficulties  through  the  blackboard  paradigm.  OPIS  uses  characteristics  of 
the  current  solution  constraints  to  focus  the  attention  of  the  top-level  manager  of  the 
knowledge  sources  (KSs).  On  each  control  cycle  OPIS  combines  constraint  propagation 
with  consistency  maintenance. 

OPIS  uses  dynamic  decompo.'^ition  to  manage  the  complexity  problem.  This  requires 
heuristic  guidance  about  the  decomposition  and  where  to  direct  the  search  effort.  Two 
decomposition  types  used  in  OPIS  are  resource-based  and  order-based.  As  implied  by 
the  naming,  a  resource-based  decomposition  emphasizes  a  specific  resource  and  stress- 
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es  optimal  resolution  of  conflicts  over  that  resource.  In  contrast,  order-based  decom¬ 
position  emphasizes  a  particular  customer's  order  and  is  directed  by  the  optimal 
resolution  of  conflicts  to  produce  that  order. 

The  architecture  of  OPIS  is  suitable  for  all  forms  of  scheduling.  In  this  framework 
knowledge  sources  are  either  analysis  or  scheduling  knowledge  sources.  Schedule  KSs 
post  updates  for  adjusting  the  incrementally  maintained  representation  of  the  current, 
incomplete  schedule.  Updates  are  one  of  these  three  kinds;  (1)  incomplete  hypotheses 
(processes  that  have  yet  to  be  schedule'^),  (2)  elementary'  conflicts,  and  (3)  opportunities 
available  due  to  loosening  of  the  constraints. 

Analysis  KSs  are  used  within  the  top-level  control  cycle,  which  consists  of  event  se¬ 
lection,  event  analysis,  action  section,  and  action  execution.  Event  selections  are  based 
on  aggregating  primitive  events  and  establishing  priorities.  Event  analysis  examines 
capacity  and  conflict  constraints.  Capacity  constraints  are  used  to  estimate  anticipated 
contention  by  identifying  possible  bottlenecks.  Conflict  constraints  are  measured  by 
duration  of  the  conflict,  the  number  of  order  involved,  projected  lateness,  etc. 

The  two  primary  schedulers  of  OPIS  are  the  order  scheduler  and  the  resource  scheduler. 
The  order  scheduler  uses  beam  search  on  resource  assignment.  The  resource  scheduler 
is  an  iterative  dispatch-based  methodology.  Other  schedulers  that  have  been  tried  sim¬ 
ply  push  scheduled  execution  times  forward  (the  right  shift  scheduler)  or  exchange  the 
remaining  part  of  an  incomplete  plan  with  another  order’s  schedule  in  order  to  optimize 
the  two  plans  (demand  swap). 


3.2  Deadline-Coupled  Real-Time  Planning 


Real-time  planning  is  an  important  class  of  planning  problems  t.hat  contain  a  mixture  of 
symbolic  and  numeric  computations.  In  real-time  planning  it  is  necessary  to  spend  time 
not  only  planning  and  acting  but  also  reasoning  about  the  time  remaining  befoie  reach¬ 
ing  a  deadline.  This  latter  reasoning  is  traditionally  done  through  meta-reasonmg,  and 
its  timing  requirements  are  frequently  ignored.  In  [Kraus  1990]  an  accounting  is  made 
for  all  planning  and  acting  times.  This  is  done  using  step-logics,  in  which  the  funda¬ 
mental  unit  (a  step)  represents  the  time  taken  to  draw  a  single  inference.  At  the  ith  step, 
the  inference  engine  uses  observations  of  the  external  environment  as  input,  together 
with  all  non-contradictory  deductions  drawn  in  the  previous  i-1  steps. 
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Until  a  plan  is  ftilly  developed,  the  reasoner  works  on  the  cuiTeni  partial  plan,  which 
consists  of  an  ordered  list  of  triples,  each  of  which  is  (1)  a  set  of  associated  precondi¬ 
tions,  (2)  an  action,  and  (3)  a  set  of  following  results.  This  representation  is  able  to 
distinguish  between  actual  facts  and  the  facts  within  the  context  of  the  partial  plan;  i.e.. 
the  difference  between  a  robot  being  actually  located  at  position  P  and  the  planning 
precondition  that  the  robot  is  located  there. 

This  approach  also  allows  execution  of  the  partially  developed  plan  to  begin  before  a 
complete  plan  has  been  constructed.  While  the  action  occurs,  planning  may  continue. 
During  acting  and  planning,  the  system  maintains  a  Working  Estimate  of  Time  (WET) 
remaining  before  the  deadline.  Funhermore,  each  action  has  a  timing  interv-al  associ¬ 
ated  with  it.  During  this  interval,  the  predicate  (action)  is  asserted  to  be  true. 


4.  Planning  Examples  in  CLP(9l) 

4,1  Departmental  Scheduling  of  Courses 

The  first  scheduling  example  is  a  program  for  constructing  a  depanment’s  schedule  of 
course  offerings  for  the  next  term.  The  program  is  somewhat  restricted,  but  does  con¬ 
tain  many  of  the  features  needed  to  assess  the  potential  of  CLP(9t)  as  an  adequate 
representation  of  scheduling  problems. 

The  entire  program  is  divided  into  three  parts:  (1)  the  rules  of  inference.  (2)  constraints 
that  must  be  met,  and  (3)  instructors'  preferences,  which  will  be  satisfied  if  possible. 
The  inference  rules  initially  construct  a  schedule  that  meets  all  constraints  and  then 
attempts  to  accommodate  each  of  the  preferences. 

It  is  assumed  that  there  is  only  one  section  for  each  of  eighteen  courses  and  that  a  par¬ 
ticular  instructor  has  already  been  assigned  to  each  course.  These  facts  are  given  in  this 
form:  course(560,jelks).  The  department  controls  only  three  classrooms  (rooml. 
room2,  room3).  All  classes  begin  on  the  hour  and  last  fifty  minutes  The  earliest  class 
begins  at  8  a.m.,  and  the  latest  class  begins  at  2  p.m.  Thus,  there  are  twenty-one 
time/room  combinations  available  for  the  eighteen  courses.  No  allowance  is  made  for 
Monday-Wednesday-Friday  or  any  other  form  of  scheduling.  One  can  assume  that  each 
class  needs  its  room  every  day. 
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Two  types  of  constraints  are  made  on  the  schedule.  First,  no  instructor  can  teach  two 
courses  at  the  same  hour.  Secondly,  certain  pairs  of  courses  cannot  be  scheduled  at  the 
same  hour  since  such  an  alignment  will  prevent  depanmenial  majors  from  taking  both 
classes.  The  second  type  of  constraint  is  represented  by  facts  like  this; 
time_clash(400,520). 

TTie  initial  scheduling  recursively  invokes  the  predicate  "scheduIe/3",  which  selects  the 
next_candidate.  This  is  an  order  scheduler  in  which  the  order  is  dictated  by  the  hour, 
with  a  numerical  sub-ordering  of  the  three  rooms.  This  scheduler  can  be  modified  to  a 
resource  scheduler.  For  example,  one  can  impose  a  priority  measure  based  on  the 
course  number  or  on  the  particular  instructor.  .A.n  alternative  is  to  construct  the  con¬ 
straint  graph  of  the  "time_clash"  predicate  and  let  the  resource  scheduler  select  the  most 
heavily  constrained  course/instructor  first.  In  the  program  of  Appendix  .-\.  the 
course/instructor  selected  for  the  current  candidate  time/room  by  the  order  scheduler  is 
based  on  a  top-to-bottom  scan  of  the  '’course"  facts. 

For  a  large  scale  system,  it  may  be  more  appropriate  to  use  a  priority  queue  of  "course  ” 
facts.  Standard  backtracking  methods  are  available  for  constructing  alternative 
schedules.  In  the  Appendix  A  example,  a  lime/room  is  left  vacant  if  no  compatible 
course/instructor  is  available.  These  unscheduled  exceptions  are  reponed  to  the  user 
when  the  Initial  Schedule  is  presented. 

After  the  initial  schedule  is  completed,  adjustments  are  attempted  for  each  instructor's 
preferences  that  are  not  satisfied  in  the  Initial  Schedule.  These  preferences  may  take 
many  forms.  An  instructor  may  wish  to  avoid  certain  hours;  e.g..  no  class  after  noon. 
She  may  wish  to  avoid  a  particular  room  assignment.  He  may  wish  to  avoid  teaching 
back-to-back  classes.  For  illustration,  the  AppendLx  ,A  program  uses  only  preferences 
for  avoiding  particular  class  hours. 

Many  priority  orders  of  the  preferences  are  available;  viz.,  by  seniority,  by  multiplicity, 
by  professorial  rank,  etc.  Each  unachieved  preference  is  tried  only  once  in  the  follow¬ 
ing  fashion.  If  (1)  there  is  an  available  time/room  in  the  current  schedule,  (2)  no 
constraints  will  become  invalid  ,  and  (3)  no  other  achieved  preference  is  negated,  then 
the  course/instmctor  is  moved  to  the  new  time/room,  and  the  old  time/room  becomes 
free.  If  any  one  of  these  conditions  is  not  met,  the  unfulfilled  preference  in  noted  (in  the 
output),  but  it  is  not  reconsidered. 


4.2  Where  Is  the  Zebra? 


The  second  example  was  originally  available  on  network  services  and  solved  by  Daniel 
Ligett.  As  stated,  one  is  given  five  house  in  a  line.  Each  house  has  one  of  five  colors 
(green,  red,  blue,  yellow,  ivory),  shelters  one  of  five  nationalities  (English,  Spanish. 
Ukrainian,  Japanese,  Norwegian),  has  one  of  five  pets  (dog.  horse,  zebra,  fox,  snails), 
stocks  one  of  five  liquids  (water,  tea,  coffee,  orange  juice,  milk),  and  breathes  one  of 
five  cigarettes  (Kools,  Chesterfields,  Parliaments,  Winstons,  Lucky  Strikes). 

With  no  constraint,  there  are  5x5x5x5x5  different  combinations  that  can  occupy  the 
first  house,  4x4x4x4x4  combinations  for  the  second  house.  3x3x3x3x3  combinations 
for  the  third  house,  2x2x2x2x2  combinations  for  the  fourth  house,  and  one  remaining 
combination  for  the  fifth  house.  Thus,  there  are  3125x1024x243x32  (about  24  billion) 
total  possible  combinations. 

Fourteen  constraints  are  given  to  limit  the  search  space  in  locating  the  house  of  the 
zebra.  All  the  restrictions  are  shown  in  the  complete  program  which  is  given  in  .Ap¬ 
pendix  B.  Eight  of  the  constraints  pair  two  of  the  five  attributes:  e.g..  the  Japanese 
smokes  Parliaments.  Three  of  the  constraints  use  the  "next_to'’  predicate;  e.g.,  the 
house  with  chesterfields  is  next  to  the  house  of  the  fox.  One  restriction  is  even  more 
specific,  using  the  predicate  ”is_to_the_right_of  ’  rather  than  "next_to."  One  constraint 
describes  the  first  house  (counting  from  left  to  right),  and  one  describes  the  middle 
house. 

With  the  fourteen  constraints,  one  finds  only  one  solution:  the  zebra  resides  in  the  fifth 
house.  If  the  constraint  that  links  the  Japanese  with  Parliaments  is  dropped,  one  finds 
nine  solutions,  and  the  zebra  can  be  found  in  any  of  the  houses  except  the  second  house. 
If  the  constraint  that  connects  snails  and  Winstons  is  also  dropped,  one  finds  sixty-four 
solutions.  If  the  constraint  between  the  green  house  and  coffee  is  dropped,  one  has  344 
solutions.  The  case  of  one  solution  and  the  case  of  nine  solutions  are  given  in  Appendix 

T3 

JL>* 

In  general,  one  cannot  expect  always  to  find  ten  solutions  by  dropping  one  of  the  four¬ 
teen  constraints.  For  example,  if  the  only  missing  restraint  is  the  one  that  requires  the 
house  of  the  Chesterfields  to  be  next  to  the  house  with  the  fox.  then  one  find  just  two 
solutions,  and  the  zebra  is  found  in  either  the  first  or  the  fifth  house. 
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The  original  solution  has  been  modified  by  adding  the  predicate  "member  restncuon.^  ' 
in  order  to  guarantee  no  duplications  of  any  of  the  five  attributes.  This  is  implied  in  the 
statement  of  the  problem  but  did  not  actually  happen  except  in  the  case  of  the  unique 
solution. 

A  variation  of  the  modified  program  also  appears  in  Appendix  B,  In  this  second  form 
the  five  values  of  each  of  the  five  attributes  is  coded  numerically.  For  example,  red- 1 . 
blue=2,  green=3,  yeIlow=4,  ivory=5  for  the  color  attribute.  Also  the  uniqueness  con¬ 
straints,  which  are  invoked  with  the  '’member_restrictions"  predicate  are  replaced  by 
two  predicates  which  use  numeric  constraints.  The  first  predicate. 
"meraber_resrrictions,"  checks  that  each  value  assumed  by  each  of  the  five  attributes 
lies  in  the  closed  interval  from  one  to  five.  The  second  predicate,  ’'column_restrictions. ' 
checks  that  no  value  is  duplicated.  If  each  numeric  value  is  assumed  only  once,  then 
their  sum  of  all  five  houses  for  each  attribute  must  be  fifteen. 


5,  Conclusions 

In  this  report  a  discussion  of  constraint  logic  programming  has  been  presented.  This 
included  (1)  features  of  the  CLP(9l)  language  that  extend  logic  programming  to  linear 
equations  and  inequalities  over  the  reals,  (2)  the  architecture  of  CLP(S^),  and  (3)  appli¬ 
cation  of  CLP(5K)  to  financial  analysis,  electrical  engineering,  and  truss  structures.  Two 
planning/  scheduling  systems  were  also  examined  to  show'  that  such  systems  could 
contain  features  that  would  be  appropriately  represented  with  CLP(3I).  Finally,  two 
small  examples  were  constructed,  a  departmental  course  scheduler  and  a  reformulation 
of  "where’s  the  zebra"  problem.  The  scheduler  illustrates  CLP(9I)'s  expressive  power 
in  planning/scheduling  problems  and  in  general  problem-solving  which  involves  sym¬ 
bolic  constraints.  The  "zebra"  puzzle  shows  the  versatility  of  CLP(SI)  with  representing 
symbolic  and  numeric  constraints. 

Future  work  in  this  area  will  be  directed  to  examining  how  the  OPIS  schedulers  and  the 
deadline -directed  real-time  scheduler  can  be  expressed  in  CLP(9l).  Hopefully,  these 
examinations  will  provide  an  impetus  for  using  CLP(9I)  in  other  areas  of  planning  and 
scheduling. 
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y*****:4e**********  +  =te*******************  +  **  +  :<:**^*:*:*******  **=»=♦**  ********/ 

/*  Scheduling  18  courses  (with  instructors)  into  three  classrooms.  */ 

/*  between  8  a.m.  and  2  p.m.  *! 

I*  Facts  are  in  files:  prog/c_and_f2.clpr,  prog/pref2.clpr  */ 

^**3tt*s|e#s(e*^*3)C3(c**  +  ***#stt5|e****5t«****=t:***********************  ******  ***=('***/ 


^*********************************************** ******* ********* 


/*  - RULES - 

/*  These  rules 

/*  go,schedule,present,time_room,  constraints_valid, 

/*  double,  overload,  and  course_clash 

/*  are  used  to  construct  the  initial  schedule,  which  honors  departmental 
/*  requirements  that  certain  pairs  of  courses  cannot  be  taught  at  the 
/*  same  hour  since  this  would  prevent  majors  from  taking  more  than  one 
/*  of  the  courses. 


*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 


^♦Jle***************************************  ******* 


:4e*******************/ 


go  schedule  ( [],  7,  room3). 
schedule(  L,  Told,  Rold) 

next_candidate(Told,Rold,Tnew,Rnew,L,Lnew), 
schedule(  Lnew,  Tnew,  Rnew). 

schedule(  L,  Told,  Rold) 

reverse(L,Lrev), 

exhibit(Lrev). 

next_candidate(Told,Rold,Tnew,Rnew,L,Lnew) 

time_room(Told,  Rold,  Tnew,  Rnew), 
course  (N,  P), 

constraints_valid(  [N,P,Tnew,Rnew],  L), 
append(  [  [NJ*,Tnew,Rnew]  ],  L,  Lnew). 

next_candidate(Told,Rold,Tnew,Rnew,L,L) 

time_room(Told,  Rold,  Tnew,  Rnew). 
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exhibit(L) 


tell(’prog/scheclule.txt’),nl,nl, 

sf: ♦  :(e +  3(c sj: :4c ♦  3|e ^ 

write(’InitiaJ  Schedule  for  CSE  Classes’),  ni, 

present(L), 

adjust_pref(L,  Ladj), 

nl,nl,write(’With  the  exceptions  listed  above,  ’), 
write(’all  preferences  have  been  met.’),nl,nl,nl, 
write('**********’*'***’*'*****************^’*'**’ ) 

write(’Adjusted  Schedule  for  CSE  Classes’),nl, 
preseni(Ladj). 


time_room(14,  roomS,  Tnew,  Rnew) 

fail. 

time_room(Told,  room3,  Tnew,  rooml) 

not(Told  =  14),  Tnew  =  Told  + 1 . 
tiine_room(Told,  rooml.  Told,  room2), 
time_room(Told,  room2,  Told,  roomS). 


constraints_valid(  [N,P,T,R],  L )  :- 

not(  member(  [N,P,_,J,  L ) ), 
not(  doubIe(  [N,P,T,R],  L  ) ), 
not(  overload(  [N,P,T,R],  L  ) ). 

double(  [N,P,T,R],  [  [N1,P1,T1,R1]  1  Tail  ] ) T=T1,  P=P1. 
double(  [N,P,T,R],  [  H  1  Tail  ] ) double!  [N,P,T,R],  Tail ). 

overload!  [N,P,T,R],  [  [N1,P1,T1,R1]  I  Tail] ) 

T  =  Tl,  course_clash!NJ^l). 

overload!  [N,P,T,R],  [  H  i  Tail] ) overload!  [N,P,T,R],  Tail), 
overload!  [N,P,T,R],  [] ) fail. 

course_clash!A,B) time_clash!A,B). 
course_cIash!A,B) time_clash!B,A). 
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present(  [  [N,P,T,R]  1  Tail  ] ) 


write(N),  write(’  ’), 
write(T),  write(’  ’), 
write(R),  wriie(’ 
write(P),  nl, 
present(Tail). 


^*sie****He5|e:t!****=<«*******5(e**5|«*5(e*****5te**  ********♦*♦  +  ♦***♦****:««♦**  ********/ 


/*  Rules  for  Adjusting  the  Initial  Schedule  to  */ 

/*  account  for  Instructors’  Time  Preferences  *I 

/*  TVL  is  the  list  of  Time  Violated  preferences  *I 


y**:*:**:t:%**:(e4:5|e*!te********=fs******************  *************  **************/ 


adjust_pref(L,Ladj) times_violate(L,[],TVL), 

adjust_tiraes(L,Ladj,TVL). 


times_violate(Lall,Lin,Lout) 

time_noz(P,T), 
not(member(  L»PiT,J,  Lin)), 
member!  L.P.T,_].  Lall), 
append!  I  L.P»T,J  ],  Lin,  Lapp), 
times„violate!Lall,Lapp,Lout). 

times_violate!Lall,Lin,Lout) 

reverse!Lin,Lout). 

adjust__times!L,Ladj,TVL) 

car!TVL,L,P,T,J), 
member!  [N,P,T,R],  L), 
reschedule!  [N,P,T,R],  L,  Lnew), 
cdr!TVL,Y),  TVLnew  =  Y, 
adjust_times!Lnew  ,Ladj  ,TVLne  w ) . 

adJust_times!L,L,G). 
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reschedule ( [N,P,T,R],  L,  Lnew) 

find_open_times(L,Lavail), 
n:ember(L._,Topen,Ropen],  Lavail), 
reraove(  [N,P,T,R],  L,  Lrem), 
check;_new_piace(N,P,T,RTopen,Ropen,Lrem,Lnew). 


find_open_times(L,Lavail) 

aIl_fimes([],7,room3,Las), 

setdiff(Las,L,Lavail). 


all_times(L,Told,Rold,Las) 

time_room(Told,Rold,Tnew,Rnew), 
append(  [  L,_,Tnew,Rnew]  ],L,Lnew), 
all_times(Lnew,Tnew,Rnew,Las). 

all_times(L,Told,Rold,Lrev) 

reverse(L,Lrev). 

check_new_place(N,P,T,R,Topen,Ropen,L,Lnew) 

constraints_valid(  [N,P,Topen,Ropen],  L), 
not(new_conflict(  [N,P,Topen,Ropen])), 
append(  [  [N,P,Topen,Ropen]  ],L.Lnew), !. 
check_new_piace(N,P,T,R,Topen,Ropen,L,Lnew) 

nl,write(’»»>  ’)» 

write(’Unable  to  accommodate  the  preference  of  ’), 
write(P),write(’  for  no  class  at  ’),write(T), 
wriie(’.’),nl, 

append(  [  [N,P,T,R]  ],L,Lnew). 

new_conflict(  [N,P,T,R]) 

time_noz(Pno,Tno),  P=:Pno,  T=Tno. 
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/*  STANDARD  USEFUL  RULES  7 

car([[N,P,T,R]lTail],[N,P,T,R]). 

cdr([[N,P,T,R]lTail],Tail). 

append([],L,L). 

append(  [XIY],  L2,  [X1L3] )  :- 

append(Y,L2,L3). 

reverse([].[]). 

reverse([HIT],L) ;-  reverse(T,Trev),  append(Trev,[H],L). 


member(  [N,P,T,R],  [  [N,P,T,R]  I J  ). 
member(  [N,P,T,R],  [  H  1  Trail  ] )  :- 

member(  [N,P,T,R],  Trail). 

remove(  [N,P,Tdl],  [  [N,P,T,R]  I  Tail  ],  Tail) !. 
removeC  [N,P,T,R],  [  H  i  Tail],  [H  1  Ldel])  :- 

remove(  [N,P,T,R],  Tail,  Ldel). 


setdiff([],B,n), 

setdiff(  [  L,_,T,R]  IA2],B,C)  :- 

meuiber(  L,_,T,R],B),  setdiff(A2,B,C). 
setdifF([  L,_,T,R]IA2],B,[  L,_XR]  IC2]) 

not(meraber(L,_XR],B)),  setdiff(A2,B,C2). 
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j*  These  are  the  course  numbers  and  the  assigned  instructor. 


course(350,deng). 

course(3 24 ,kabuid) . 

course(360,elalemam). 

course(400,faisal). 

course(520,han). 

course(624,baker). 

course(660,kabuki). 

course(200,amm). 

course(220,baker). 

course(405,govel). 

course(505,govel). 

course(524,han). 

course(530,ijssel). 

course(533,amm). 

course(560,jelks). 

course(561,jeiks), 

course(618,deng). 

course(630,faisal). 

/*  These  are  the  restrictions  on  scheduling  courses  at  the  same  hour.  */ 

/*  Ultimately  in  the  rules  each  combination  is  forbidden.  */ 

time_clash(200,350). 
time_clash(220,350). 
time_clash(324,350). 
time_clash(324,360). 
time_clash(360,400). 
time_clash(400,520). 
time_clash(405,530). 
time_clash(505,533). 
time_clash(520,560). 
time_clash(524,56 1 ). 
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time_clash(6 1 8.624). 
time_clash(624,660). 
time_clash(630,660). 


/*  These  are  instructors’  preferences,  constraints  which  can  be  relaxed. 
/*  time_noz(A,B)  means  instructor  A  teaches  at  B  o’clock. 

/*  room_noz(A,B)  means  instructor  teaches  in  room  B. 

J*  adj_noz(A)  means  instructor  has  adjacent  or  back-to-back  classes. 
!*  Ultimately  the  rules  will  deny  all  three  of  these  preferences. 

time_noz(amin,8). 
time_noz(amin,  14). 
time_noz(kobuki,  12). 
time_noz(faisaJ,10). 
tLme_noz(govel,  10). 
time_noz(govcl,ll). 
time_noz(govel,  1 2). 
time_noz(ijssel,8). 
time_noz(jelks,8). 
time_noz(jelks,9). 

room_noz(  deng  ,r  3 ) . 
room_noz(ijssel,r  1 ). 

adj_noz(baker). 
adj_noz(deng). 
adj_noz(e  lalemain) . 
adj_noz(han). 


APPENDIX  B 
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THE  ORIGINAL  VERSION  OF  ’’WHERE’S  THE  ZEBRA  ” 

/*  Each  house  is  represented  by  a  structure  of  the  form: 
house(Colour,  Nationality,  Pet,  Drink,  Cigarette) 

Puzzle  is  currently  missing 

next_to(house(_, chesterfields),  house(_._»fox,_,_).  Houses)  */ 

puzzle 

houses(Houses), 

member(house(red,  english,  _,  _),  Houses), 
member(house(_,  Spanish,  dog,  _),  Houses), 
member(house(_,  ukranian,  tea,  _),  Houses), 
member(house(_,_,snails,_,winstons).  Houses), 
member(house(green,_,_, coffee Houses), 
member(house(_, Japanese, parliaments).  Houses), 
right_of(house(green,_._,_,J,  house(ivory Houses), 
meraber(house(yellow,  _,  kools).  Houses), 

Houses  =  L.  houseC,  _>  niilk,  _),  _,J, 

Houses  =  [houseC.  norwegian,  _)IJ, 

next_to(house(_, kools),  houseC._.liorse,__,_),  Houses), 

meraber(house(_,  _,  _,  orange  Juice,  lucky  _strikes).  Houses), 

next_to(house(_,norwegian,_,_,_),  house(bIue,_, _,_._),  Houses), 

member_restrictions(Houses), 

print_houses(Houses), 

nl  nl  ^i'ite(’***************’<'************************’),nl,fail. 
houses([ 

house(_,  -X 

house(_. 

houseC,  _), 

house(_.  _). 

houseC.  _>  _) 

]). 
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member_restrictions(Houses) 

member(house(green,_, Houses), 
member(house(red, Houses), 
member(house(blue,_,_,_,_).  Houses), 
member(house(yellow, Houses), 
member(house(ivory Houses), 
inember(house(_,english,_,_,_),  Houses), 
member(house(_,spanish,_,_,_),  Houses), 
member(house(_,ukxanian,_,_,_).  Houses), 
member(house(_,japanese,_,_,_).  Houses), 
member(house(_,norwegian,_,_,_).  Houses), 
member(house(_,_, horse ,_,_),  Houses) , 
ineinber(house(_,_,ciog,_,_),  Houses), 
member(house(_,_, zebra, _,_),  Houses), 
member(house(_,_,fox,_,_),  Houses), 
member(house(_,_, snails, _,_),  Houses), 
member(house(_,_,_,orangeJuice,_),  Houses), 
member(house(_,_,_,niilk,_),  Houses), 
member(house(_,_,_,coffee,_).  Houses), 
member(house(_,_,_,tea,_).  Houses), 
member(house(_,_,_,  water, _),  Houses), 
member(house(_,_,_,„,kools).  Houses), 
member(house(_,_,_,_, chesterfields).  Houses), 
member(house(_,_,_,_.,lucky_sn’ikes).  Houses), 
member(house(_,_,_,_,parliaments).  Houses), 
niember(house(_,„,_,_,winstons).  Houses). 


right_of(A,  B,  [B,  A  I  J). 

right_of(A,  B,  [X  i  Y]) right_of(A,  B,  Y). 

next_to(A,  B,  [A,  B  I  J). 
next_to(A,  B,  [B,  A  i  J). 
next_to(A,  B,  [X  I  Y]) next_to(A,  B,  Y). 
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member(X,  [XIY]). 

member(X,  [AIB]) member(X,  B). 

print_houses([AlB3) tell(’prog/zout.txt’), 
print(A),nl, 
print_houses(B). 
print_houses(0). 


The  Output  When  All  14  Restriciton  Are  Included 


house(yeIlow,  norwegian,  fox,  water,  kools) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(red,  english,  snails,  milk,  winstons) 
house(ivory,  Spanish,  dog,  orange_juice,  lucky_strikes) 
house(green,  Japanese,  zebra,  coffee,  parliaments) 
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The  Output  When  Only  13  Restriciton  Are  Included 
Missing:  house(_japanese,_,_, parliaments) 

house(yellow,  norwegian,  fox,  water,  kools) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(red,  englisli  snails,  milk,  winstons) 
house(ivory,  Spanish,  dog,  orangejuice,  lucky_strikes) 
house(green,  Japanese,  zebra,  coffee,  parliaments) 


:1c  :tc  :1c  :tc  :1c  :1c  :1c  :1c  :1c  :1c  :1c  :1c  :1c  :1c  :1c  :1c  :|c  :1c  :1c 

house(yeIlow,  norwegian,  fox,  water,  kools) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(red,  english,  zebra,  milk,  parliaments) 
house(ivory,  Spanish,  dog,  orangejuice,  lucky_strikes) 
house(green,  Japanese,  snails,  coffee,  winstons) 


house(yellow,  norwegian,  zebra,  water,  kools) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(red,  english,  fox,  milk,  parliaments) 
house(ivory,  Spanish,  dog,  orangejuice,  lucky _strikes) 
house(green,  Japanese,  snails,  coffee,  winstons) 

^  ak  •ly  at*  a^  a,^  -jr 

^a  a^  a^  7^  ^fC  a^  ^C 

houseCyellow,  norwegian,  fox,  water,  kools) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(red,  english,  snails,  milk,  winstons) 
house(ivory,  Japanese,  zebra,  orangejuice.  lucky _strikes) 
house(green,  Spanish,  dog,  coffee,  parliaments) 


♦  :<s  **  =(!*=♦:*★*****:<:  5(! ^  :(c  :|e  :(e 

house(yellow,  norwegian,  zebra,  water,  kools) 
house(blue,  ukranian,  horse,  tea,  parliaments) 
house(red,  english,  snails,  milk,  winstons) 
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house(ivory,  Japanese,  fox,  orange  juice,  lucky  _srrikes) 
house(green,  Spanish,  dog,  coffee,  chesterfields) 


3fc  :<c  4:  :(e  :4: 4c  :1c  :1c  :4e  :tc  :t:  ^  ^  ***  ^  ********  4c  *:<(**  ^ 

house(yellow,  norwegian,  zebra,  water,  koois) 
house(bIue,  Japanese,  horse,  orange  Juice,  iucky_strikes) 
house(red,  english,  snails,  milk,  winslons) 
house(ivory,  ukranian,  fox,  tea,  parliaments) 
house(green,  Spanish,  dog,  coffee,  chesterfields) 


house(yellow,  norwegian,  fox,  water,  koois) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(ivory,  Spanish,  dog,  milk,  parliaments) 
houseCgreen,  Japanese,  snails,  coffee,  winstons) 
house(red,  english,  zebra,  orange  Juice,  lucky_s  trikes) 

**************************************** 

houseCyellow,  norwegian,  fox,  water,  koois) 
house(blue,  ukranian,  horse,  tea,  chesterfields) 
house(ivory,  Japanese,  snails,  milk,  winstons) 
houseCgreen,  Spanish,  dog,  coffee,  parliaments) 
house(red,  english,  zebra,  orange  Juice,  lucky  _strikes) 


**************************************** 

houseCyellow,  norwegian,  zebra,  water,  koois) 
houseCblue,  ukranian,  horse,  tea,  parliaments) 
houseCivory,  Japanese,  snails,  milk,  winstons) 
houseCgreen,  Spanish,  dog,  coffee,  chesterfields) 

houseCred,  english,  fox,  orange  Juice,  lucky_strikes) 
**************************************** 
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THE  NUMERICAL  VERSION  OF  '’WHERE’S  THE  ZEBRA” 


/*  Each  house  is  represented  symbolically  by  a  structure  of  the  form;  *l 

I*  house(Color,  Nationality,  Pet,  Drink,  Cigarette)  */ 


1* 

TRANSLATIONS 

COLORS 

NATIONALITffiS 

PETS 

DRINKS 

SMOKES 

l=red 

l=english 

l=dog 

1=  water 

l=kools 

2=blue 

2=spanish 

2=:horse 

2=iea 

2=chesteifields 

3=green 

3=:ukranian 

3=zebra 

3=coffee 

3=parliamenis 

4=yellow 

4=japanese 

4=fox 

4=orangeJ 

4=winstons 

5=ivory 

5=norwegian 

S-snails 

5i=milk 

5=lucky_  strike 

/^Missing  constraint  of  Japanese  and  parliaments:  member(house(_,4,_,_,3),Houses)*/ 
puzzle 

houses(Houses), 

member(house(l, _),  Houses), 
member(house(_,  2, 1 ,  _),  Houses), 

member(house(_,  3,  2,  _),  Houses), 

member(house(_,_,5,_,4),  Houses), 
member(house( 3  ,_,_,3 ,_),  Houses), 
right_of(house(3, house{5, Houses), 
member(house(4,  _,  1),  Houses), 

Houses  =  L.  houseC,  5,  _),  _,  J, 

Houses  =  [houseCj  5,  _)U, 

next_to(house(_, 2),  house(_,_,4 ,_,_),  Houses), 
next_to(house(_,_,_,„,l),  house(_,_,2,_,_).  Houses), 
member(house(_,  „,  _,  4,5),  Houses), 
next_to(house(_, 5 house(2,_,_,_,_).  Houses), 
member_restrictions(Houses), 
column_restrictions(Houses), 
print_houses(Houses), 
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house(_,  J, 

house(_,  _)^ 
house(_,  _)« 
house(_,  _). 
houseL,  _) 


]). 


meinber_restrictions(Houses) 

member(house(CJS[,P,D,S),  Houses), 
1  <=  C,  C  <=  5, 

1  <=  N,  N  <=  5, 

1  <=  P,P<=5, 

1  <=  D,  D  <=5, 

1  <=  S,  S  <=  5,  fail. 
member_restriction(Houses). 


column_restrictions(Houses) 

Houses  =  [  house(Cl,Nl,Pl,Dl,Sl), 
house(C2,N2,P2,D2,S2), 
house(C3.N3,P3,D3,S3), 
house(C4,N4,P4,D4,S4), 
house(C5,N5,P5,D5,S5)  }, 
Cl+C2+C3+C4+C5=15, 

N  1+N2+N3+N4+N5= 1 5, 
Pl+P2+P3+P4+P5=15, 
Dl+D2+D3+D4+D5=15, 
Sl+S2+S3+S4+S5=15. 

right_of(A,  B,  [B,  A  I  J. 

right_of(A,  B,  [X  I  Y  ] right_of(A,  B,  Y). 

next_to(A,  B,  [A,  B  I  J). 
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next_to(A,  B,  [B,  A  i  J). 

next(to(A,  B,  [X  I  Y]) next_to(A,  B,  Y). 


raember(X,  [XIY]). 

member(X,  [AIB]) meraber(X,  B). 

print_houses([AiB]) 

teil(’prog/zout.txt’), 
print(A),  nl, 
print_houses(B). 

priiit_houses([]). 


"US  GOVfBNMfNT  PBINTINf,  Of  net 
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MISSION 

OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  te  nology  transition  in  support  of  Air 

3 

Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C  I  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/ maintainability  and  testability. 


